Low temperature structural model of hep solid C70 
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We report intermolecular potential energy calculation for solid C70 and determine the optimum static 
orientations of molecules at low temperature which are consistent with the monoclinic structural model 
proposed by us in an earlier report [Solid State Commun., 105, 247 (1997)]. This model indicates that the 
Cs-axis of the molecule is tilted by an angle ~ 18° from the monoclinic b-axis in contrast with the molecular 
orientation proposed by Verheijen et al.[J. Chem. Phys., 166, 287 (1992)] where the Cs-axis is parallel to 
the monoclinic b-axis. In this calculation we have incorporated the effective bond charge Coulomb potential 
together with the Lennard- Jones potential between the molecule at the origin of the monoclinic unit cell 
and its six nearest neighbours, three above and three below. The minimum energy configuration for the 
molecular orientations turns out to be at 6 = 18°, 
orientations. 

PACS nos.: 61.48. +c, 61. 50. Ah, 61.50.Lt 



= 8° and if) = 5° where 9, <j> and tp define the molecular 



The orientational ordering transformation in solid Cgo 
is well understood!! Potential energy calculations using 
only the Lennard- Jones potential derive the low. tempera- 
ture structure of solid Ceo to be orthorhombiccl, in sharp 
contradiction with the experimentally observed simple 
cubic structure. Lu et al.B could successfully derive the 
simple cubic structure for the low temperature phase of 
solid Ceo by adding a Coulomb potential term to the 
Lennard-Jones potential. This Coulomb interaction be- 
tween neighbouring molecules arises from the effective 
charges at the single bonds and the double bonds of 
the molecules in the ordered phase at low temperature. 
In the case of solid C70, the elongated molecular shape 
causes much greater diversity in the structural phase 
transformation sequence as well as in the orientational 
ordering of molecules. One of the first experimental stud- 
ies of structural phase transformation, on cooling, in solid 
C70 is by Verheijen et alfl and they reported the following 
sequence of phase transformations 

fee # rh <- hcp-1 (c/a « 1.63) # hcp-2 (c/a w 1.82) 
^ monoclinic 

where the fee is the high temperature phase and the 
monoclinic is the low temperature phase. The lattice pa- 
rameters reported for the monoclinic unit cell are, a = 
c = 19.96 A , b = 18.51 A and (3 « 120° with space 
group P112i/m. They proposed a structural model for 
this monoclinic structure in which the orientation of the 
long axis (C5) of the static molecules is parallel to the 
monoclinic b-axis (i.e., the original hep c-axis). Subse- 
quently, they have also carried out the potential energy 
calculation using the atom-atom Lennard-Jones poten- 
tial and the corresponding result,was in agreement with 
their model. Later, Que et al.El pointed out that the 
unique angle f3 for a monoclinic structure with space 
group P112i/m cannot be 120°, but should be less than 
120° which is not consistent with the molecular orien- 
tations proposed by Verheijen et al. They have carried 
out a detailed group theoretical analysis of librations in 
solid C70 and proposed an orthorhombic structure for 



which the two space groups Pbnrn or Pnma have the same 
Raman modes. Nelissen et aLQ have carried out a har- 
monic lattice dynamical calculation using the Lennard- 
Jones potential which proposed that the structure of solid 
C70 should be monoclinic at low temperature with space 
group P2i/m. This structure requires two independent 
angles which they defined as 01 and <j>2 and pointed out 
that for their structure <f)l + <fi2 should be 36°. Later, 
Agterberg et al.13 have pointed out that the calculation of 
Nelissen et al. can also admit of an orthorhombic struc- 
ture with space group Phnm. In another detailed study 
using the anisotropic interaction potential between two 
C70 molecules in terms of Cayley-Klein parameters, they 
have confirmed that the structure Wuld be orthorhom- 
bic with space group Pbcm ot Pbnnia. The experimental 
results of van Smaalen et al.El also conclude in favour of 
an orthorhombic structure with space group Pbnm of the 
solid C70 at low temperature. It is to be noted that this 
structure has resulted by a two step transformation from 
the ideal h.c.p. phase at high temperature. Oh et al.t2l 
have reported detailed potential energy calculations us- 
ing Girifalco-Lennard- Jones potential for different molec- 
ular orientations depicting various possible structures of 
the solid and the corresponding lattice parameters. They 
have compared their results with experimental reports by 
othersE2l. They reported two possible transformation se- 
quences on cooling, viz., 

fee ^ ortho(I,II) ^ rh <— hep I ^ mono 

or, fee <— hep II ^ hep (III, IV) ^ mono 
where the fee is the high temperature phase and the 
monoclinic is the low temperature phase. Here hep I - 
IV correspond to the c/a ratios of the unit cells which 
are 1.84, 1.61, 1.50 and 1.53 respectively. The ortho I 
and II correspond to the lattice parameters, a = 16. 2A , 
b = c = 14.0A and a = 16.0A , b = 14.4A , c = 
14. 8A respectively. The monoclinic structures in both 
cases are the same and correspond to the Cs-axis of ori- 
entation being parallel to the b-axis of the unit cell. This 
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is similar to the molecular orientation in the monoclinic 
structure predicted by Verheijen et alia, but the lattice 
parameters in .the two reports are significantly different. 
Pickholz et al.Ell, while using atom-atom Lennard-Jones 
interaction in their calculation, obtain the low temper- 
ature structure to be hep with c/a ~ 1.77, instead of 
a monoclinic structure. The molecular orientations in 
their model-are, however, parallel to the unit cell c-axis. 
Sprik et alJl3 have carried out constant pressure molec- 
ular dynamics calculations incorporating Lennard-Jones 
potential together with the bond charge Coulomb po- 
tential in a phase transformation study of fee solid C70. 
Their result indicates a two step transformation, viz., fee 
^ rh monoclinic which is in excellent agreement with 
the experimental results reported by Christides et al.ES. 
In the low temperature monoclinic phase predicted by 
these authors, the molecules in alternate layers of the 
structure are oriented parallel to the (111) direction and 
to the (110) direction respectively. We have reported in 
our earlier papenlil that the ideal hep (c/a ~ 1.63) solid 
C70 undergoes a single step transformation, on cooling, 
to a monoclinic phase which is in contrast with the two 
step transformation reported by Verheijen et alJa. We 
conjectured that the single step transformation implies 
simultaneous freezing of the molecular rotations around 
their long axes (Cs-axes) and the short axes (C^-axes), 
as compared to sequential freezing, which results in a two 
step transformation. The former is likely to lead to a dif- 
ferent monoclinic structure as compared to the latterQ'El. 
We had proposed a structural model, consistent with the 
lattice parameters for the monoclinic structure measured 
from our x-ray diffractogram (XRD) at 100 K, wherein 
the molecular Cs-axis has a tilt of about 18° from the 
monoclinic b m -axis. Our recent experimental studies oa. 
the kinetics of phase transformation in hep solid C70C3 
show that, except for very slow cooling rates, a two step 
transformation invariably results. Such a two-step trans- 
formation may be similar to that reported by Verheijen 
et alo and Christides et al.E3. The low temperature mon- 
oclinic structures and the C70 molecular orientations re- 
ported by these authors must, therefore, correspond to 
some local potential energy minima of the system raised 
above the global minimum by misoricntation stresses. 

In this paper, we report the results of the potential en- 
ergy calculation pertaining to the monoclinic structure 
obtained by us from our 100 K XRD. Since the cool- 
ing rate was very slow (typically, 0.0033 K /min.) the 
monoclinic structure determined from our experiment is 
expected to be close to the equilibrium structure at low 
temperature. In this calculation, we have considered the 
interactions between the seven molecules, one at the ori- 
gin (0,0,0) of the monoclinic unit cell and the remaining 
six at the six nearest neighbour positions, three above 
and three below the basal plane. The orientation of the 
molecules has been defined using the angles, 8, cj) and ip, 
as explained later. The final orientation of the molecules, 
as obtained by locating the minimum of the potential en- 
ergy has 9 = 18°, = 8° and V = 5°. 



The structure of our C70 sample at room temperature 
is ideal hcplij and the XRD pattern of this structure has 
been analysed in detail in an earlier workllZl. The hep 
lattice constants are, a/ t = b^ = 10.53A , Ch — 17.26A 
and 7/j = 120°. The 100 K pattern fits well to a mon- 
oclinic structure with lattice parameters, a m = 10.99A , 
b m = 16.16A Jj m = 9.85A and f3 m — 107.75° and space 
group P2\ /11a It is clear that the following correspon- 
dence of lattice parameters is appropriate in going from 
the hep to the monoclinic: Ch shrinks and forms b m , b^ 
increases and forms a m , and a^ shrinks to form c m . This 
correspondence is shown in Fig.l. The hexagonal angle 
7h can be seen to reduce to the monoclinic unique angle 
P m by a shear along the direction shown by the arrows in 
Fig. 1(b). It can be seen that the above correspondence 
indicates a tilt of the Cs-axis of C70 molecules by an angle 
of approximately 18° from the b m -axis of the monoclinic 
unit cell (see Fig. 1(b)). 

In this calculation we have used the Lennard-Jones po- 
tential together with the bond charge Coulomb potential. 
Following Lu et al.0 in the ordered state of the solid there 
is less charge density distribution around the C-C single 
bonds (long bonds) compared to that around the C-C 
double bonds (short bonds). Such a charge density dif- 
ference on the surface of the molecule leads to short range 
intermolecular Coulomb interaction which decays as 1/ 
R 5 at long distances in the case of C70. In this calcula- 
tion we have used the same notation for the interaction 
sites as was done by Sprik et al.E3 in their calculation for 
the fee solid C70. The Lennard-Jones centres are taken 
to be the carbon atoms sites (C sites: ace — 3.4A ) as 
well as the mid-points of the electron-rich C-C double 
bonds (D sites: (Jdd = 3.6A , whereas <jcd — 3.5A ) 
and the magnitude, of the interaction has been taken as, 
e = 2.964 meVBEI 

The effective charge around the singleJaonds are as- 
sumed to be localised around the C sitestij. Though in 
C70 molecules, there are five types of C sites, for the sake 
of simplicity all C sites can be treated to be of the same 
type as in the case of Ceo- If a C site carries an effective 
charge qc , the corresponding value at the D sites will be 
q_D = - 2qc- The values of these effective charges are 
fractions of tht; electronic charge, e.g., qc? = 0.175e and 
qo = - 0.35eE3. The bonds around the equatorial region 
of the molecule have a length intermediate between the 
length of a single bond and a double bond and, therefore, 
are termed as intermediate bonds (or, I bonds). For these 
bonds, <7j = 3.5A and qj = 0.5qj> Since the polar cap of 
the molecule consists of both long bonds and short bonds, 
the Coulomb interaction is stronger at the two opposite 
poles of a molecule in the ordered phase. Therefore the 
bond charge interaction between two molecules in adja- 
cent planes is dominated by the polar caps and can be 
termed " polar cap" interaction. The coordinates of the 
70 carbon atoms in the molecule weie generated from 
data available in published literature £9. The total inter- 
action potential between any two C70H molecules is given 
by 
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where the subscript x indicates interaction sites, e.g., C, 
D, or, I, and rij,2j are the coordinates of the atoms, 
t>im,2n coordinates of bond centres and q m ,n are effec- 
tive bond charges (e.g., qc, qz? and qj). 

Using the above interaction potentials between the 
seven nearest neighbour molecules for the lattice param- 
eters measured from our 100 K x-ray diffractogram, we 
optimise the molecular orientations with respect to the 
unit cell parameters by minimising the total energy. The 
molecular orientations are defined using the polar coor- 
dinates 9 and (f>, where 9 is the polar angle, i.e., the angle 
between the b m -axis and the molecular long axis and <p is 
the azimuthal angle. Further minimisation of the poten- 
tial energy has been performed by varying the angle (V>) 
which defines the body rotation (spin) of the molecules 
around their Cs-axes of the six molecules above and be- 
low the basal plane relative to the molecule at the origin 
of the monoclinic unit cell. 

We calculate the potential energy per cluster for each 
configuration of the seven molecules defined by the val- 
ues of 6 between 0° and 90°, <f> between 0° and 180°, 
and ifj between 0° and 72°. In Fig. 2, we show the vari- 
ation of potential energy with respect to 9, <f> and ip. 
Clearly, for = 8° and ip = 5° the potential energy is 
minimum at 9 = 18°. This result is consistent with the 
model we proposed for the molecular orientation in the 
monoclinic unit cell experimentally measured by ust-X 
The minimum potential energy per cluster obtained for 
our model is -3.923 eV which is significant^ lower than 
other calculations reported in literatureallll It is to be 
mentioned here that there are two other minima with 
respect to 9 at 54° and 88° (not shown in the figure). 
Though these three minima have energy values compara- 
ble to each other, the energy barriers between them are 
nearly 0.04 eV which is too large to spontaneously flip 
the molecule from one minimum to the other at low tem- 
peratures (100 K). These multiple minima indicate that 
the system may have degenerate states with respect to 9, 
or, that they may even correspond to different structures 
which can be obtained by applying some external driv- 
ing force, e.g., pressure etc. The degeneracy may arise 
as a consequence of the repulsive interactipas between 
the C-C bonds as indicated by Sprik et aB. in their 
extended model. Since the present calculation is based 
on the static orientations of molecules (molecular static 
calculation) the study of the structural transformation 
or the change of molecular orientations with tempera- 
ture cannot be carried out by this method. A detailed 
molecular dynamics calculation is needed for this pur- 
pose which is in progress. Nevertheless, we argue that 
the equilibrium low temperature structure will be con- 
sistent with the molecular Cs-axis tilted away from the 



unit cell long axis, because of the polar cap interaction, 
as proposed in the present model and supported by the 
potential energy calculation, rather than being oriented 
parallel to itu'cl. We recall from the introductory sectios. 
that the molecular dynamics calculation of Sprik et al.c3 
leads to two competing orientations for the molecules. A 
pictorial view of the optimised static orientation of the 
seven molecules for our monoclinic lattice parameters is 
shown in Fig. 3. This picture is generated using a public 
domain software MOLDRAW. 

In conclusion, our results of potential energy calcula- 
tion using the Lennard - Jones potential together with 
the bond charge Coulomb potential support the-,struc- 
tural model proposed by us in our earlier papenU. The 
calculation shows that the minimum potential energy 
state appears with a tilt of the Cs-axis of the molecule 
with respect to the monoclinic long axis (b m -axis) by an 
angle of 18°. Further minimisation of the interaction en- 
ergy between the neighbouring molecules has been car- 
ried out using an angle V which defines the rotational 
configuration of the molecule around its own Cs-axis. 
The molecular orientations obtained from our calcula- 
tion have an implication of significant polar cap interac- 
tion between the nearest neighbour molecules in adjacent 
planes. The bond charge interaction is thus seen to play 
a key role in determining the equilibrium low tempera- 
ture structure of solid C70. Our calculations show also 
that the structural model proposed by us corresponds to 
the most stable low temperature structure of solid C70 
reported so far. 
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FIG. 1. (a) High temperature hep structural model; C70 
molecules are freely rotating mimicking spheres, (b) Mono- 
clinic unit cell model proposed for the lattice parameters ob- 
tained from XRD pattern at 100 K (Ref. 10). 

FIG. 2. Variation of the potential energy in eV/ cluster (a) 
with 6, (b) with (f>, and (c) with ip 

FIG. 3. Final orientational configurations of the 7 
molecules of the unit cell arising from the potential energy 
calculation; a m , b m , and c m represent the monoclinic cell pa- 
rameters. 
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